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摘要 : 太阳 图 像 中 包含 明显 的 太阳 活动 区 域 ， 以 及 无 现象 的 太阳 宁静 区 。 从 太阳 图 像 中 
识别 有 效 的 太阳 活动 区 域 ， 是 图 像 处 理 技术 在 天 文 研究 中 的 典型 应 用 。 得 益 于 美国 太阳 动力 
学 天 文 台 的 日 球 物理 学 事件 知识 库 提供 的 实时 太阳 观测 数据 ， 提 出 了 一 种 基于 日 球 物理 学 事 
件 知识 库 的 太阳 活动 识别 方法 。 此 方法 获取 6 种 太阳 活动 的 信息 (发 生 时 间 、 位 置 、 区 域 面 
积 ) ， 建 立 对 应 时 间 的 全 日 面 图 像 的 尺度 变换 模型 。 结 合 位 置 与 区 域 面积 信息 ， 对 不 同 太阳 
活动 进行 梯度 阅 值 分 割 ， 边 界 识别 方法 被 用 来 定位 和 识别 太阳 活动 的 区 域 。 然 后 ， 对 太阳 图 
像 特征 参数 相关 性 的 研究 ， 得 到 每 种 太阳 活动 的 最 住 特征 参数 组 合 。 方 法 实现 了 对 太阳 活动 
的 精确 定位 和 有 效 识 别 ， 为 后 续 工作 的 开展 提供 了 便利 。 此 外 ， 对 不 同 太阳 活动 区 域 提取 特 
定 组 合 的 特征 ， 可 以 为 基于 内 容 的 图 像 检索 系统 建立 精简 的 图 像 特征 集 提供 一 种 可 行 办 法 。 
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基于 内 容 的 图 像 检 索 技 术 产 生 于 20 世纪 90 年 代 ， 它 的 出 现 解决 了 大 规模 数字 图 像 检 索 的 难题 。 
随后 ， 一 大 批 研究 性 或 商用 的 基于 内 容 的 图 像 检 索 系 统 ( Content-Based Image Retrieval, CBIR ) 应运 而 
生 。 区 别 于 原 有 系统 中 对 图 像 进行 人 工 标注 的 做 法 ， 基 于 内 容 的 检索 技术 自动 提取 每 幅 图 像 的 视觉 内 
容 特征 作为 索引 ， 如 色彩 、 纹 理 、 形 状 等 。 而 针对 大 规模 太阳 活动 的 基于 内 容 的 图 像 检 索 系 统 也 已 出 
现 ， 如 太阳 动力 学 天 文 台 的 基于 内 容 的 图 像 检 索 系 统 (SDO Content-Based Image Retrieval System) 。 

虽然 现 有 的 面向 天 文大 数据 的 基于 内 容 的 图 像 系 统 能 够 满足 在 海量 天 文 图 像 库 中 检索 目标 图 像 的 
基本 需求 ， 但 是 依然 面临 检索 耗 时 、 准 确 率 不 高 的 问题 。 出 现 这 些 问题 的 原因 如 下 : (1) 原始 采集 的 
太阳 图 像 尺 寸 大 ， 对 应 所 提取 的 特征 点 明显 增多 ， 造 成 耗 时 增加 ; (2) 太阳 图 像 中 呈现 混淆 难 辨 的 太 
阳 现 象 ， 所 提取 的 特征 没有 针对 性 ， 造 成 匹配 结果 不 理想 。 本 文 的 研究 对 以 上 问题 提供 了 解决 思路 。 

本 文 的 目标 在 于 为 建立 太阳 活动 自动 识别 和 检测 系统 探讨 可 行 的 办 法 ,该 系统 的 功能 是 自动 识别 
一 张 给 定 的 全 日 面 太阳 图 像 中 所 包含 的 太阳 活动 。 建 立 这 样 一 个 系统 所 面临 的 难题 : (1) 不 同 太阳 活 
动 的 物理 特征 差异 很 小 ， 且 同 种 太阳 活动 包含 多 种 物理 形态 ,使 得 太阳 活动 的 判定 不 明确 ; (2) 目前 
针对 多 种 太阳 活动 识别 的 技术 方法 还 不 成 熟 ; (3) 特征 库 匹 配 时 间 开 销 大 ， 检 索 效率 低 。 

针对 问题 (1) ， 本 文采 用 日 球 物理 学 事件 知识 库 提供 的 实时 太阳 观测 数据 ， 数 据 由 太阳 动力 学 天 
文 台 的 特征 检测 小 组 (Feature Finding Team, FFT) 模块 自 动 获 取 ， 避免 了 人 为 干预 ， 数据 包括 太阳 实 
时 活动 的 发 生 时 间 、 在 太阳 表面 的 发 生 位 置 以 及 活动 类 型 ， 针 对 问题 (2) ， 本 文采 用 网 格 结构 划分 图 
像 ， 使 用 梯度 闷 值 法 分 割 太 阳 活 动 的 目标 区 域 ， 使 用 边界 识别 方法 最 小 外 接 和 矩 形 ( Minimum Bounding 
Rectangle, MBR) 和 八 连通 链 码 获取 目标 区 域 的 轮廓 特征 数据 ;针对 问题 (3) ， 本 文 研究 了 6 种 图 像 参 
数 (主要 是 纹理 特征 ) 所 属 特 性 ， 探 讨 图 像 参 数 与 太阳 活动 的 相关 性 ， 找 出 每 种 太阳 活动 对 应 的 最 优 
参数 或 参数 组 合 ， 并 依次 对 包含 6 类 太阳 活动 的 目标 区 域 图 像 集 进行 特征 提取 ， 用 于 建立 特征 库 。 
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1 相关 工作 


太阳 动力 学 天 文 台 '" 是 美国 国家 航空 航天 局 一 个 为 期 10 年 的 太阳 观测 任务 。 太 阳 动 力学 天 文 台 
的 科学 目标 是 在 小 尺度 的 时 间 和 空间 下 以 多 波段 研究 太阳 大 气 层 ， 从 而 了 解 太阳 对 地 球 和 近 地 球 太空 
区 域 的 影响 。 大 气 成 像 组 件 ( Atmospheric Imaging Assembly, AIA) 用 于 拍摄 高 时 间 与 高 空间 分 辩 率 的 完 
整 太阳 影像 。 事 件 检 测 系 统 ( Event Detection System, EDS) 利用 特征 检测 小 组 自动 产生 特征 数据 ， 并 传 
送 到 日 球 物理 学 事件 知识 库 ( Heliophysics Event Knowledgebase，HEK)。 文 [2] 利 用 分 层 迭 代 法 以 及 最 
优化 方法 对 日 球 物理 学 事件 知识 库 提供 的 太阳 观测 数据 进行 太阳 活动 间 的 相关 性 研究 ， 最 终 求解 网 络 
流 。 目 前 ， 该 方法 仅 对 太阳 活动 区 和 上 日冕 洞 两 种 太阳 活动 的 追踪 效果 较 好 。 文 [3] 结 合 日 球 物 理学 事 
件 知 识 库 的 太阳 观测 数据 ， 对 方形 活动 区 域 提取 10 种 特征 参数 ， 构 造 用 于 太阳 活动 区 域 检 索 的 特征 
描述 符 ， 并 研究 了 对 应 太阳 活动 的 最 佳 相 似 匹配 方法 。 结 果 显 示 ， 对 于 除 太 阳 活 动 区 、 耀 班 和 S 型 结 
构 3 种 太阳 活动 外 ， 该 方法 还 不 能 获取 满意 的 检索 准确 率 。 
针对 单一 太阳 活动 现象 的 自动 识别 成 为 非常 流行 的 研究 课题 。 文 [4] 利 用 神经 网 络 贝 叶 斯 
= (Bayes) 推理 和 形状 关系 自动 识别 包括 太阳 黑子 、 耀 班 和 日 冕 物质 抛射 。 文 [5] 使 用 局 部 阔 值 和 区 域 
> 生长 法 在 全 日 面 图 像 中 识别 太阳 活动 区 。 文 [6] 训练 了 神经 网 络 识 别 太阳 图 像 的 隐藏 背景 中 的 暗 条 。 
文 [7] 提 出 了 一 种 自动 的 暗 条 检测 方法 ， 这 种 方法 超越 了 典型 的 暗 条 检测 ， 它 确定 暗 条 的 疹 和 方向 
fa, 发现 暗 条 手 性 。 文 [8] 提 出 了 一 种 自动 识别 次 斑 算 法 ， 这 种 算法 在 噪声 和 背景 有 亮 线 影响 的 情况 
下 效果 很 好 。 对 于 日 冕 物质 抛射 的 识别 ， 文 [9] 提 出 了 一 种 基于 数据 立方 体 的 算法 ， 并 利用 小 波 变换 
分 析 识别 日 震 。 目 前 的 方法 还 仅 限于 识别 单一 太阳 活动 或 物理 特征 相近 的 几 种 太阳 活动 ， 适 用 性 不 
- 佳 ， 无 法 满足 海量 天 文 图 像 处 理 的 需求 。 
m 虽然 ， 目 前 的 识别 算法 结合 了 太阳 活动 的 物理 特征 ， 并 综合 运用 了 图 像 处 理 的 各 种 技术 ， 因 此 对 
m 单一 太阳 活动 识别 取得 了 比较 好 的 结果 ， 但 是 还 没有 通行 的 方法 针对 全 日 面 图 像 上 的 多 种 太阳 活动 进 
N 行 识 别 。 随 着 天 文 图 像 数 量 的 激增 ， 现 有 的 太阳 活动 识别 方法 在 时 效 性 和 识别 精度 方面 不 能 满足 多 波 
~ 段 太阳 图 像 的 识别 任务 。 而 日 球 物理 学 事件 知识 库 实时 提供 太阳 活动 的 精准 数据 源 ， 这 些 数据 源 记录 
了 太阳 活动 物理 特征 的 重要 参数 ， 可 以 为 相关 研究 工作 提供 有 力 的 参考 (0 。 现 有 的 方法 中 结合 日 球 
物理 学 事件 知识 库 与 图 像 处 理 方法 用 于 太阳 活动 定位 和 识别 的 相关 研究 还 不 多 。 本 文 介绍 一 种 基于 日 
球 物理 学 事件 知识 库 的 太阳 活动 定位 和 识别 的 方法 ， 以 实现 对 全 日 面 图 像 太 阳 活 动 进行 自动 识别 。 
本 文 从 日 球 物理 学 事件 知识 库 获 取 6 类 太阳 活动 的 实时 观测 数据 ， 包 括 活动 区 ( Active Region, AR), 
© H ja] (Coronal Hole, CH) , 7 FRX (Emerging Flux, EF) , WP (Flare, FL), F (Sunspot, SS), S 
型 结构 (Sigmoid, SG) ， 相 比 其 它 类 型 的 太阳 活动 ， 以 上 6 类 太阳 活动 发 生 频 繁 ， 且 持续 时 间 长 。 
太阳 动力 学 天 文 台 位 于 太空 中 的 成 像 系统 大 气 成 像 组 件 获 取 不 同 波段 的 紫外 线 光谱 成 像 ， 拍 摄 的 
太阳 图 像 不 受 大 气 透 明度 和 大 气 扰动 等 因素 的 影响 ， 但 是 太阳 全 日 面 上 的 亮度 在 紫外 线 光 谱 通道 边缘 
变 亮 ， 出 现 临 边 增 亮 现象 ( Limb Brightening) 。 本 文 首 先 对 太阳 图 像 进行 标准 化 处 理 ， 去 除 太阳 边缘 
亮 圈 ， 避 免 非 日 面 区 域 对 后 续 工 作 的 影响 。 然 后 ， 从 日 球 物理 学 事件 知识 库 采 集 365 张 全 日 面 图 像 的 
所 有 太阳 活动 的 相关 信息 ， 包 括 太阳 活动 的 发 生 时 间 、 发 生 位 置 、 覆 盖 面 积 等 ， 为 了 将 这 些 信 息 应 用 
于 标准 化 后 的 全 日 面 图 像 ， 需 要 建立 对 应 时 间 的 全 日 面 图 像 的 尺度 变换 模型 ， 产 生 的 新 位 置信 息 和 面 
职 信息 将 应 用 于 本 文 的 图 像 分 割 方法 中 。 
采用 梯度 阔 值 法 ， 结 合 不 同 太阳 活动 的 物理 特性 与 尺度 变化 后 的 新 数据 ， 对 不 同 太 阳 活 动 分 区 域 
进行 分 割 ， 将 分 割 后 的 太阳 活动 目标 区 域 以 两 种 边界 识别 方法 标注 显示 ， 记 录 目 标 区 域 的 轮廓 信息 。 
为 了 建立 用 于 基于 内 容 的 图 像 检 索 系统 的 精简 图 像 特 征集 ， 本 文 研究 了 太阳 图 像 6 种 特征 的 相关 
性 。 对 2015 全 年 的 全 日 面 太 阳 图 像 提 取 6 种 特征 参数 ， 通 过 分 析 6 种 特征 相关 性 ， 找 出 相关 性 最 强 
的 图 像 特征 组 合 后 就 能 实现 降低 特征 参数 的 维度 。 结 合 目标 区 域 的 轮廓 信息 ， 提 取 目 标 区 域内 的 特征 
组 合 ， 建 立 特征 集 。 
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本 文 提出 一 种 基于 日 球 物理 学 事件 知识 库 的 太阳 活动 识别 方法 ,该 方法 结合 日 球 物理 学 事件 
知识 库 提供 的 太阳 活动 观测 数据 ， 利 用 多 种 图 像 处 理 方法 ， 实 现 了 对 尺度 可 变换 的 全 日 面 图 像 中 多 种 
太阳 活动 的 准确 定位 与 识别 。 实 验 结果 表明 ， 该 方法 在 太阳 活动 识别 准确 度 及 时 间 开 销 方面 取得 满 
意 的 效果 。 


2 太阳 活动 识别 方法 


2.1 Hough 变换 检测 圆心 和 半径 

Hough 变换 是 模式 识别 领域 对 二 值 图 像 进行 圆 检测 的 有 效 方法 ， 本 文 提出 一 种 利用 圆 斜 率 特性 的 
改进 Hough 变换 方法 ， 用 于 检测 太阳 圆心 和 半径 以 及 去 除 边缘 亮 圈 。 362: I8 B JS EE ae s 间 域 变 
换 到 参数 空间 ， 用 大 多 数 边界 点 满足 的 某 种 参数 形式 描述 图 像 中 的 区 域 边界 。 通 过 设置 累加 器 进行 累 
加 ， 求 得 峰值 对 应 的 点 就 是 所 需要 的 信息 。 
2.1.1 圆 的 几何 特性 分 析 

圆周 上 存在 4 个 点 PCs, y). POR Yuin) > POR Yes POGQ S yw) 分 别 对 应 圆 的 左 端 、 上 
端 、 右 端 、 下 端的 极 值 点 ， 除 这 4 个 极 值 坐 标点 外 ， 对 于 圆周 上 的 任意 一 点 POS, y), ， 过 点 P 了 存在 一 
条 垂直 方向 的 线段 六 与 一 条 水 平方 向 的 线段 请 ， 广 、 旋 与 圆 的 交点 分 别 是 4、 有 8，AP4B 是 过 点 P 的 
内 接 直角 三 角形 ， 根 据 圆 的 几何 特性 : 圆 的 内 接 直角 三 角形 的 斜 边 过 圆心 ， 因 此 ，48 线段 的 中 点 就 
是 圆心 点 。 限 定 P 的 移动 范围 ， 构 造 4 种 圆 内 接 直 角 三 角形 ( 如 图 1): 

(a) A {P| xi, LE < Aa V LY L Yndi CDA {Plann «x < Xas Yan LY < HHS 

(c) 为 {Plx, «x «x, Yu LY SYY; (DOH {Plea «x «x, y, «Y « Yno 


WAAL 


Al 4 种 圆 内 接 直角 三 角形 


Fig. 1 Four kinds of circle inscribed right-triangle 


ix B(x,, yi), 则 PB 的 中 点 坐标 为 xu=(x + x1) /2, yo=(y + y1)/26 
2.1.2 Hough 变换 检测 圆心 和 半径 
(1) 利 用 索 贝尔 (SOBEL) 算 子 从 原始 图 像 中 检测 出 边缘 信息 ， 并 转化 边界 为 单一 像素 的 边缘 图 ， 
然后 将 图 像 二 值 化 ( 即 设 定 边缘 像素 值 为 1， 其 他 像素 值 为 0) ， 将 边缘 像素 坐标 存 入 数组 V(x, y)。 
(2) 设 置 4 个 累加 需 数组 C(x, y) (初始 值 为 0)， 按 照 2.1.1 节 方 法 ,逐个 判断 V(x, y) P P 
属 情况 ， 计 算 中 点 坐标 (xu，yo) ， 并 将 对 应 情况 的 累加 器 数组 C(x, y) 加 1。 累 加 器 数组 C 中 的 最 大 
值 对 应 的 坐标 即 为 圆心 (p, 9) 。 
(3) 将 圆心 坐标 (p, gq) 和 V(x, y) 代 入 圆 方程 ， 利 用 Hough 变换 对 参数 R 进行 统计 ， 当 RR 的 统计 
EKF 80%( 即 V(x,y) 中 80% 的 点 落 在 圆周 上 )，R 即 为 圆 的 半径 。 
(4) 根 据 圆心 和 半径 ， 将 图 像 超出 半径 的 区 域 设置 为 背景 ， 实 现 去 除 边缘 亮 圈 的 目的 。 
2.2 图 像 建 模 方法 
2.2.1 尺度 变换 
对 365 张 原始 全 日 面 图 像 进 行 标准 化 处 理 后 ， 图 像 的 大 小 尺度 发 生变 化 ， 除 太阳 活动 的 发 生 时 间 
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7 不 变 外 ， 太 阳 活 动 的 发 生 人 位置、 覆盖 面积 随 图 像 尺度 大 小 的 变化 而 变化 。 设 原始 采集 信息 中 太阳 活 
动 的 发 生 位 置 为 (x_ pre, y. pre) , Eust f47J Area _ pre， 图 像 中 太阳 半径 为 R_ pre, 经 2.1 市 标准 化 
处 理 后 , 设 太阳 活动 的 新 位 置 为 (x_ new, y. new), AWIN Area_ new， 图 像 中 太阳 半径 为 R_ new, 
变换 过 程 如 下 : 

如 图 2(a) ,假设 太阳 活动 的 发 生 位 置 位 于 太阳 边缘 点 ， 那 么 存在 以 下 关系 : 


ls (R. new/R., pre) (1) 
y. new = y. pre (R_new/R_pre) | 
如 图 2(b) ,假设 太阳 活动 为 阴影 部 分 ， 其 面积 5 的 计算 如 下 : 
S = [ 2mrRdo = | 2mR?eosbdg = 2m? | cos0d0 = 2mR' (1 ~ sind) = 27RH . (2) 
0 6 0 


R_pre 


(x_pre, y_ pre) 


X new, y_ new) 


. pre 


一， 
E 
Ve 

bw 


(a) 
图 2 尺度 变换 的 几何 模型 


Fig.2 Geometric model for scale transformation 


由 (2) 式 推导 太阳 活动 的 覆盖 面积 : 


Area. pre = 27R_ preH 
(3) 


Area. new = 2mR. newh 


H : x 
Aks faee = ， 进 而 推导 出 如 下 关系 : 


Area, pre 2 


=K". (4) 


Area _ new 


由 2. 1 节 得 到 尺度 变化 后 的 新 图 ， 建 立 用 于 图 像 分 割 的 直角 坐标 系 (图 像 的 首 个 像素 点 作为 坐标 
系 原点 ) ， 坐 标 系 变化 后 ， 对 太阳 活动 的 位 置 坐标 (x, y) 作 再 次 调整 


x. new - R, new * x, new 


y_new = R new - y new (5) 
2.2.2 网 格 划分 
网 格 划 分 是 常规 的 图 像 处 理 方法 ， 对 标准 化 处 理 后 的 图 像 用 4 x 4 像素 的 网 格 进行 划分 ,考虑 到 
低 维 的 数据 结构 要 比 高 纬 的 数据 结构 在 运算 效率 方面 表现 更 优 ， 本 文采 用 之 前 研究 中 的 处 理 方法 ， 对 
网 格 数据 进行 降 维 操作 和 存储 '""] 。 


Chinay nA 4eE RAT 
ChinaX ive ERBTI 


2 期 李 卫 疆 等 ， 基 于 日 球 物理 学 事件 知识 库 的 太阳 活动 识别 255 


2.3 边界 识别 方法 
最 小 外 接 和 矩形 和 八 连通 链 码 是 图 像 分 类 和 识别 研究 常用 的 方法 ， 被 用 来 提取 目标 区 域 的 几何 特 
JEU, AS SCH) SP BEET EB TE Ad 23-96 BE BIER LRL, ESSE [HEU] 6 种 太阳 活动 区 域 。 特 
别 地 ， 最 小 外 接 矩 形 信息 存储 在 POINT (xi, Yi, X2, Yar Wax Y3, Was Yes Xis Yı) m , 而 链 码 信息 采用 
POLYGON (x yo X4, Y2, 05, +085 X,, Vn» O) ETM, FLA (x, y, wo) 表 示 网 格 首 像素 的 坐标 及 
方向 。 
2.4 BEBE BEARS 
本 文采 用 梯度 阔 值 分 割 法 ， 结 合 6 种 太阳 活动 的 物理 特性 ， 在 满足 判定 条 件 的 情况 下 ,分 配 不 同 
BS BUE, WHE 2.2.1 节 新 坐标 系 中 太阳 活动 位 置 坐 标 (x_new, y_ new) ， 将 该 点 的 灰 度 值 作 为 初始 靖 值 ， 
对 图 像 的 网 格 单元 进行 初步 分 割 ， 利 用 网 格 单元 的 最 小 外 接 和 矩形 修正 闷 值 。 梯 度 国 值 分 割 法 过 程 如 下 : 
输入 : 一 张 经 过 Hough 变换 及 尺度 变换 处 理 后 的 图 像 IMG 
输出 : 一 张 所 有 太阳 活动 区 域 被 标注 后 的 图 像 
步骤 : 1; Convert IMG to a gray image 
2: IMG is divided into M rows and N columns of cells 
3: IMG—Block ( pixel, row, pixel, col, cell, row, cell_col) 


— Block ( pixel, row, pixel, col, index) 


4: Block( pixel, row, pixel, col, index) —img 


eo 5: Set img(x new, y, new) to the initial threshold T, 
J 6: if ~ exist MBR( img) 
一 set threshold to T,( T, « Tj) 
l ue 7: else if exist MBR(img) A Area( MBR) / Area_new>threshold_area 
set threshold to T,(T,> T,) 
一 8 : else exist MBR(img) A Area( MBR) / Area_new<threshold_area 


set threshold to T,( T, < T,« T,) 
9: count the number; [ Block( pixel, row, pixel col, index) »T] (Te[T,, Ti, T,, T,]) 


= 10: if T > threshold_number 
T 11; set Block( pixel. row, pixel col, index) to 255 
= 12; else set Block( pixel row, pixel col, index) to 0 
c 13; count the density of connected square grids 
14; if density < threshold_density 
15; remove the grid area 
16; else end 


经 过 Hough 变换 、 尺 度 变 换 等 处 理 后 ， 根 据 太阳 活动 的 位 置信 息 ， 将 新 的 图 像 分 区 域 进行 网 格 划 
分 ， 并 依次 对 6 种 太阳 活动 进行 梯度 国 值 分 割 ， 使 用 最 小 外 接 和 矩形 和 八 连通 链 码 标注 太阳 活动 区 域 。 其 
中 ， 步 又 5 将 太阳 活动 的 位 置 点 的 灰 度 值 作为 初始 阔 值 7,， 步 又 6 设置 判定 条 件 (1): 若 初 次 分 割 后 的 
识别 区 域 不 存在 最 小 外 接 和 矩形 ， 设 定 新 闻 值 7,(7T,< T) 。 步 又 7 设置 判定 条 件 (2) : 若 初次 分 割 后 的 
识别 区 域 存 在 最 小 外 接 和 矩形 ， 并 且 Area( MBR) 49 Area _ new 面积 的 比值 大 于 经 验 闵 值 threshold _ area, 
设 定 新 闻 值 7,(7,> T) 。 步 又 8 设置 判定 条 件 (3) : 若 初 次 分 割 后 的 识别 区 域 存 在 最 小 外 接 抢 形 ， 并 
且 4rea(MBR) 与 4rea_new 面积 的 比值 小 于 经 验 靖 值 threshold_area， 设 定 新 间 值 7,(7Ti< T3< T,)o 7b 
又 6~8 根据 3 种 判定 条 件 ， EMRE BEC < Tu< 7,< T.) 。 对 于 每 种 太阳 活动 ， 以 上 3 种 判定 条 件 
可 灵活 调整 ， 以 及 分 配 不 同 的 梯度 冰 值 。 本 文中 ，threshold _ number 取经 验 阔 值 0.15, threshold _ 
density LAI EY {A 0. 5, threshold_areaC[0.3, 0.7]。 不 同 太 阳 活 动 的 物理 特性 存在 差异 ， 太 阳 活 动 
DC, FK, WEBEL S 型 结构 相 较 于 黑子 、 日 罗 洞 ， 亮 度 更 加 明显 (次 斑 是 最 显著 的 太阳 活动 ， 黑 
子 、 日 金 洞 亮度 则 很 低 )。 因 此 ， 选 取 合 适 的 立 值 对 太阳 图 像 分 割 产生 有 益 的 效果 。 
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2.5 ”特征 参数 相关 性 评价 方法 X1 图 像 特征 参数 列表 
首先 介绍 6 种 图 像 寺 征 参数。 #1 中 .L Table 1 List of Extracted Image Parameters 
为 太阳 活动 区 域 的 像素 个 数 ，z 表示 第 i 个 像 。 编号 名 称 
素 的 灰 度 值 ，p(z,) 表示 灰 度 = 在 灰 度 直方 图 。 
中 的 比值 。 
对 太阳 图 像 自动 特征 提取 方法 已 经 有 了 部 p oe 
AVES C873. xx e JE IE MZ Bh EL EB BC A 
中 筛选 表现 优异 的 参数 组 合 ， 而 涉及 各 特征 参 p3 e 
数 间 关 联 性 的 研究 则 较 少 。 为 了 研究 各 特征 参 
数 间 的 关联 性 ， 本 文采 用 如 下 方法 。 P4 偏 度 系数 dd b evite 
(1) 提 取 太 阳 图 像 活动 区 域 X( 外 接 和 矩形 
的 标注 范围 ) 的 6 种 特征 参数 ，g,(X,) 表 示 提 ps 均匀 性 U= Ý p) 
取 第 大 个 活动 区 域 的 第 i 种 参数 。 B ; 
(2) 计 算 太 阳 活 动 区域 X Pi PSM P6 和 
第 j 个 参数 的 协 方差 矩阵 ; 
1 N D, 
< 2 [g;(X,) - ui] * [q;(X,) -u] , (6) 
CO — JU. us LY Dg); N=365; D, H n 张 太阳 图 像 所 包含 的 活动 数量 。 
e (3) 计 算 各 参数 间 的 相关 矩阵 R,， 反 映 了 第 i 和 第 j 两 个 参数 间 的 相关 程度 ， 即 R, 越 大 ， 两 特征 
= 越 相关 。 
N RS, (7) 
: TOD 
3 35 验 
3.1 数据 集 


本 文 图 像 数据 集 来 源 于 太阳 动力 学 天 文 台 的 大 气 成 像 组 件 (2015 4F)®, FE 365 张 太阳 全 日 面 图 
像 ， 图 像 基本 信息 如 表 2。 
表 2 太阳 活动 的 基本 信息 


Table 2 A summary of event types 


太阳 活动 来 源 报道 次 数 采集 模块 
太阳 活动 区 193° A 20226 空间 可 能 性 聚 类 算法 模块 
H Aid 193^ A 16 590 空间 可 能 性 聚 类 算法 模块 
MEDE 131° A 21774 Y SAT eo fh o Y 
S 型 结构 131°A 8 538 S 型 结构 检测 模块 
浮现 区 日 震 与 磁 成 像 仪 10 467 浮现 区 检测 模块 
黑子 日 震 与 磁 成 像 仪 4 493 欧洲 太阳 黑子 观测 模块 


太阳 动力 学 天 文 台 通 过 仪器 捕捉 不 同 波段 的 太阳 紫外 线 ， 由 大 气 成 像 组 件 进行 光谱 成 像 ， 合 成 具 
有 明显 太阳 活动 的 图 像 ， 不 同 太 阳 活 动 被 事件 检测 系统 的 相应 模块 报道 。 例 如 ， 太 阳 黑 子 由 日 震 
与 磁 成 像 仪 (Helioseismic and Magnetic Imager, HMI) 检测 ， 这 种 现象 在 2015 年 全 年 被 EGSO_SFC 模块 
报道 了 4493 次 。 


(D http: //sdowww.|msal.com/suntoday_v2/ 
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3.2 太阳 活动 区 域 定位 与 识别 

本 文采 用 3 种 波段 紫外 线 ( 波 长 211、193、171) 著 加 后 形成 的 图 像 ， 这 些 图 像 中 太阳 活动 区 、 焰 
PE, Hæ, FAK 4 种 太阳 活动 的 物理 特征 较为 清晰 。 利 用 日 球 物理 学 事件 知识 库 提 供 的 信息 ， 结 
合 多 种 方法 ， 在 全 日 面 图 像 上 完成 了 太阳 活动 的 定位 和 识别 工作 。 最 小 外 接 和 矩形 准确 定位 了 完整 的 太 
阳 活 动 区 域 ， 八 连通 链 码 较 好 地 展示 了 太阳 活动 的 形态 特征 。 如 图 3。 


| Events | 


(b) 


图 3 (a) 来 源 于 大 气 成 像 组 件 的 全 日 面 图 像 ; (b) 太阳 活动 分 割 与 识别 
Fig.3 (a) Original image of the Sun taken on 1 January 2015 from AIA; (b) Segmentation and Recognition of solar activities 


3.3 识别 准确 度 评 估 

本 文 进一步 统计 了 2015 年 1 月 31 天 的 6 种 太阳 活动 的 分 制 比例 (分 割 比例 是 指 八 连通 链 码 识别 
区 域 面 积 与 最 小 外 接 和 矩形 定位 区 域 面积 的 比值 ) ， 用 于 评估 太阳 活动 识别 的 准确 度 。 上 具体 过 程 : 采用 
3 种 方法 (种 子 区 域 生长 法 (SRG)、 分 水 岭 分 割 法 (WaterShed ) IAE BE B] (B. 43 BITE (GTS) ) 对 太阳 活 
动 定位 区 域 进行 活动 识别 ， 然 后 ， 分 别 统计 3 种 方法 应 用 于 图 像 集 后 的 分 割 比 例 。 结 果 见 图 4 和 图 5。 
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图 4 识别 准确 度 比 较 


Fig.4 Comparison of recognition accuracy of three methods 
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图 5 时 间 开 销 比较 


T Fig. 5 Comparison of time consumption of three methods 

N 图 4 对 比 了 3 PII ARAE, ADR AMAF 186 个 太阳 活动 区 域 ， 对 离散 数据 点 的 进 
co 一 步 分 析 发 现 ， 使 用 种 子 区 域 生 长 法 得 到 的 186 组 数据 中 ， 共 有 70 组 分 割 比例 在 0.3 ~0.7 范围 内 ， 
c 其 它 116 组 数据 是 由 于 图 像 过 度 分 割 或 者 图 像 空洞 等 原因 产生 的 ， 这 些 数据 被 认为 无 法 满足 识别 要 
NL 求 ， 采 用 种 子 区域 生 长 法 的 识别 准确 率 为 37.6%。 相 应 地 ， AKI UU RIERA 48. 9% ( SE 
m 有 91 STATS BER), BABE BIA BUTI TA MEM Sim, y 65. 6% (FEA 122 组 数据 符合 要 求 ) 。 


XT A438] LC Bs EC SEES AR, GBD 2. 4 Trys obe BE pO vie ZS o 
图 5 比较 了 3 种 分 割 方法 的 时 间 开 销 ， 其 中 值得 注意 的 是 ， 太 阳 活 动 区 以 及 日 冕 洞 这 两 种 太阳 活动 

的 发 生 范围 较 大 ， 因 而 在 时 间 开 销 上 大 于 其 余 6 种 太阳 活动 。 统 计 结果 显示 ， 种 子 区 域 生长 法 的 平均 时 

间 开 销 为 0. 111s， 分 水 岭 分 割 法 和 梯度 阀 值 分 割 法 的 平均 时 间 开 销 分 别 为 0. 049 s, 0.035 s。 

3.4 特征 参数 相关 性 评价 
由 2.5 节 方 法 对 2015 年 1 月 31 天 的 太阳 图 像 进行 特征 相关 性 计算 ,将 6 种 太阳 活动 图 像 各 特征 

参数 相关 性 结果 用 热度 图 表示 ( 白色 为 1， 黑 色 为 0) ， 如 图 6。 


E 将 6 种 太阳 活动 不 同 特征 间 相 关 性 进行 可 视 化 表征 ， 得 到 以 下 结果 : 太阳 活动 区 图 像 的 可 选 特征 
o 参数 组 合 包括 (P1, PA) , (P2, P3), (P4, P5, P6); E GILL A RT XE RH CUR (P1, P5, 
P6), (P4, P5) ; 浮现 区 图 像 的 可 选 特征 参数 组 合 包 括 (P1, P3), (P2, P4, P6) ; 次 斑 图 像 的 可 选 特 


征 参 数组 合 包 括 (P1, P2), (P3, P4, P6) ; 黑子 图 像 的 可 选 特征 参数 组 合 包括 (P4，P5，P6) (P3, 
P4); S 型 结构 图 像 的 可 选 特征 参数 组 合 包括 (P1, P3, P5), (P4, P5) 。 对 太阳 活动 区 域 提取 特征 参 
数 ， 应 优先 选择 耗 时 最 少 的 组 合 。 


4 总 结 


/UN 


本 文 从 太阳 动力 学 天 文 台 的 日 球 物理 学 事件 知识 库 采 集 6 种 太阳 活动 的 基本 信息 ， 通 过 对 图 像 建 立 
尺度 变换 模型 ， 结 合 每 种 太阳 活动 的 物理 特征 对 图 像 活动 区 域 进行 梯度 闪 值 分 制 ， 使 用 最 小 外 接 矩 形 和 
弗 里 曼 ( Freeman) 码 实现 对 太阳 活动 区 域 的 定位 与 识别 。 对 太阳 图 像 各 特征 参数 的 相关 性 进行 研究 ， 
得 出 了 合适 的 参数 组 合 。 研 究 结果 说 明基 于 日 球 物理 学 事件 知识 库 的 太阳 活动 识别 方法 的 可 行 性 。 但 
仍 有 可 改进 的 地 方 ， 主 要 是 对 各 太阳 活动 物理 特征 在 灰 度 空间 的 表现 ， 需 要 加 深 了 解 ， 有 助 于 改进 梯 
度 阔 值 分 割 方法 ， 进 一 步 提高 对 太阳 活动 识别 的 准确 度 。 未 来 工作 将 对 太阳 活动 区 域 提取 特定 的 组 合 参 
数 ， 应 用 精简 数据 集 建立 基于 内 容 的 图 像 检 索 系统 ， 并 测试 在 检索 耗 时 与 检索 准确 率 等 方面 的 表现 。 
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Recognition of Solar Activities Based on Heliophysics Event Knowledgebase 
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Technology of Kunming, Kunming 650500, China, Email; aininalal @sina.cn) 


Abstract; Solar image conveys obvious active regions and peaceful regions of no activity. Recognizing 


useful solar activities from solar image is a typical application of image processing in the astronomical research 


scope. Benefiting from real-time observation data provided by Heliophysics Event Knowledgebase ( HEK ) of 


Solar Dynamics Observatory ( SDO), this paper proposes a recognition method of solar activities based on 


HEK. In this method, we gather six kinds of solar activities datas ( date, location, area) , then we design a 


scaling transform model for full-disk solar image of corresponding date. Combining location data and area data, 


we design different gradient thresholds to segment different solar activities regions. Two kinds of boundary 


identification method are used to locate and recognize solar activities regions respectively. Furthermore, we 


study the correlation of six kinds of solar image features, and give combinations of features for every solar 


activity. The method in the paper realizes the precise positioning and efficient recognition, which lies the 


foundation for further work. Besides, extracting particular combination of features from different solar activities 


regions also provides a practical approach of constructing a condensed feature set for CBIR system. 


Key words: Solar activity location; HEK; Gradient thresholds; Combination of features 


